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In this paper, we present the recent progress on minimally doubled lattice actions. In particular, 
we discuss the proposal of Creutz and its variations on an orthogonal lattice. A preliminary 
computation of the pion mass on an SU(3) background field shows the expected behaviour as 
predicted form the chiral perturbation theory. 
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1. Introduction 

Motivated by the Dirac structure of the graphene electrons in two dimensions, Creutz intro- 
duced a lattice action, which in four dimensions, describes two flavours of Dirac fermions with 
exact chiral symmetry In momentum space, this action has the form: 

3 

D{p)=iBY4i^-Y,cosp^) + i Y^yMp) , (1-1) 

where 

si (p) = sin pi + sin p2 — sin pi, — sin p4 
Slip) = sin;?! — sin p2 — sin pi, + sin p^ 
sj,{p) = sinpi — sin p2 + sin p^, — sin p^ 

and B,C are free parameters. Creutz has shown that the zeros of this operator in the Brillouin zone 
are at {p,p,p,p) and {—p,—p,—p,—p), where C = cos^. In order to get a lattice action with 
one zero at the origin, we proposed a formulation, which shifts the zeros in the Brillouin zone at 
(0,0,0,0) and {2p,2p,2p,2p) [§. 

The minimal number of zeros allowed by the Nielsen-Ninomiya theorem is exactly two 
Hence, these actions achieve a minimal doubling of fermion species on the lattice. In fact, the idea 
of minimally doubled action is not new. It was first pointed out by Karsten and later by Wilczek 
[@, |5|]. Their actions are unitary equivalent to each other. In the free case, the Wilczek action has 
the form: 

4 3 
D{p)= 'T^M sin;?^ + /74 52(1- cos ;?i) . 

It has two zeros: one at (0,0,0,0) and the other at (0,0,0, n). 



2. The Creutz Action on Orthogonal Axes 

The easiest way to get a lattice action with one zero at the origin is to put the Creutz action on 
an orthogonal lattice, where the action parameters satisfy BS = C, where S = sin p. Setting B = I 
and translating momenta, p^ = p + q^, we get: 

3 3 
D{q) = /74(4C - Cj^cos^^) + iC ^ Jk^kiq) + iSy^Ymiq^ + iS ^ JkCkiq) , 

H k=l H k=\ 

where we have denoted: 

ci {q) = (coS(7i — 1) + (cosg2 — 1) — (cosg3 — 1) — (cosg4 — 1) 
C2{q) = (cos<7i — 1) — (cosg2 — 1) — (cos 173 — 1) + (cosg4 — 1) 
C'i{q) = (cos^i — 1) — (cos ^2 — 1) + (cosg3 — 1) — (cos^4 — 1) 

Taking C = 1 / \/2 as well as 5 + C = 0, and defining: 

Si\{q) = — sin ^1 — sin ^2 — sin ^3 — sin ^4 

c/^{q) = (cos<7i — 1) + (cosg2 — 1) + (coS(73 — 1) + (cosg4 — 1) , 
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the p-translated Creutz action takes the form: 

D{i) = Y^^ytiSii{i)+Y^iyiiCti{q) , 

or in the scalar product notation, {y,x) = Yii^ii, one has: 

D{q)=i{Y,s{q)+c{q)) . 
Using the following orthogonal matrices: 

/I 1 -l-l\ 

11-1-11 , 1 

a := — b := — 

2 1-11-1' 2 

V-1 -1 -1 -l) 
and noting that, 

s = las, c = 2bc , 

where 

5 = ( sin^i, sin(72, sin(73, sin 174 c = {cosqi — \, cos(72 — 1, cos(73 — 1, cos(74 — 1 
then, the rescaled action by a factor of 2 can be written in the form: 

D{q) := i{Y,aS{q) +bc{q)) = i{a^Y}^{q) +a^bc{q)) . 



(I 1 -1 -l\ 
1-1-1 1 
1-11-1 

VI 1 1 ly 



Denoting, 



a'.= a 



/-II 1 l\ 

I- 111 

II- 11 

VI 1 1 -ly 



we get: 



D{q) = i{a''Y,s{p) + ac{q)) . 
It is easy to show that a^/ are Dirac gamma matrices :Pl and Pr 

{(fl'^7)/ii(«^7)v} = Y^ap^QaviYp^Ya} = 2Y^ap^apv =25^v . 
p,a p 

Therefore, the factor can be dropped. This way, the final expression has the form: 

D{q) = /(7,%)) + /(/,%)) 

= L iJp. sin q^+Y, iYfii^o^ q^-'^) , 

where / = a/ are again Dirac gamma matrices for the same reason as above. Noting that, 

we get another expression for the fermion action: 

D{p) = '/m + L ~ • 

This expression was elegantly derived by Creutz in terms of a linear combination of two naive 
actions plus the —2iY term, the latter cancelling exactly one naive action at the zeros of D [0]. 
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2.1 The Dirac Operator in Position Space 

In order to write down the Dirac operator in position space, we express the momentum space 
operator in terms of forward and backward propagating plane waves, e'Pf^ and e^'Pf^ : 

Dip) ='n + ^l^[{y^-ir^,) e'P^ + (/^ + ir^,) e-'P^] - liT , 

where we have added the bare fermion mass, m. Then, by making the formal substitution ip^ — > d^, 
one gets: 



D = ml + - X (/^ - iy^) A + (/^ + iy^) e 



- 2/r , 

where the shift operators, e^>^ and e^^>^ , are defined by their action on a Dirac field ^{x): 



Hence, the position space Dirac operator can be implemented using the following terms: 

• the on-site term, (ml — 2/r)i//,-. 

• the forward hopping term, ^(Y^ — iy^) + afi). 

• the backward hopping term, ^{Y^ +iy^)Y{x — afi). 

As usual, the gauge fields are introduced by requiring the hopping terms to be gauge covariant. 

3. Minimally Doubled Actions and Hypercubic Symmetry 

Shortly after our proposal, it was noted that the lattice action given above lacks the full sym- 
metry of the hypercubic group The reason is that the action picks as a special direction the main 
diagonal of the hypercube. Hence, in the presence of gauge field interaction, there is a dimension 
five operator that enters the action, namely \j/rV^Y [@]- 

As we pointed out earlier, the idea of minimally doubled action is not new. The actions of 
Karsten and Wilczek pick the time axis as a special direction. Again, the loss of the full hypercubic 
symmetry introduces extra relevant terms in the interacting case ^ ^. 

Since a minimally doubled action has necessarily two zeros, the line that joins the two ze- 
ros in the Brillouin zone defines a special direction. Hence, independently of a specific action, 
the hypercubic symmetry will be broken. In principle, there is nothing special about a hypercu- 
bic action. The authors of reference ^ propose a 'hyperdiamond' action, which generalises the 
graphene structure in five dimensions. However, the resulting action has more than two zeros. It 
is an interesting question whether it is possible to have a minimally doubled action which does not 
break the original symmetry of the action. 

4. A Preliminary Numerical Test 

So far, the present formulation has only been tested at tree level perturbation theory [Q]. It 
has been shown that the scaling violations are of the order 0{a^), as expected. Here, we present 
preliminary results on the pion mass calculations on ten 16^32 lattices generated with the SU(3) 
Wilson gauge action at /3 = 6. 
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4.1 The Quark Propagator 

The quark propagator, g, is computed as the solution of Unear system Dg = 5, where 5 is taken 
to be a point source. Since the massless operator is antihermitian, it is easy to see that D is normal, 



i.e. D*D = DD* . This property, which is shared by the Kogut-Susskind operator []12|], allows one to 
use optimal inversion algorithms. Indeed, the Conjugate Gradients algorithm on Normal Equations 
(CGNE) and the Conjugate Residual (CR) algorithm are optimal Krylov subspace based algorithms 



for staggered fermions [|10|]. Hence, the minimally doubled fermion presented here shares the same 
numerical advantages, while describing two species of fermions instead of four. For the quark 
propagator computations we have used the CGNE algorithm. 
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Figure 1. The scaling law of CGNE iterations with the bare quark mass. The critical exponent 
is computed to be z = 0.92(4). 

While optimal, this algorithm suffers from the critical slowing down, as all Krylov subspace- 
based inversion algorithms do. Thus, the the number of CGNE iterations to reach a fixed accuracy 
is expected to scale like ~ 1/ {amy with the inverse quark mass, am. Our data suggest a critical 
exponent value z = 0.92(4), which is clearly smaller than 1, the expected value for Wilson fermions 



[10] 



In the present calculations, we haven't made any effort to accelerate the inversion. The even- 
odd acceleration does not work for these type of fermions since the on-site term does not commute 
with the hopping term. 

4.2 The Pion Propagator and the Pion Mass 

Zero momentum pion propagators are computed using the vacuum expectation values of cor- 
relation functions of pion interpolating fields, ^[x)ys^f{x): 

G{t) = £ < o\w)75w{xm^)75w{m > 
= '£g{x,org{x,o) . 
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For large t the right hand side can be fitted to the ground sate ansatz: 

G{t) ~cosham;r(r/2 + l-f) , 

where periodic boundary conditions are applied and T is the lattice extension along the fourth 
diection. We computed effective masses by inverting the expression: 

G{t + \) cosham;r(r/2-f) 



G{t) cosham;r(r/2 + l-0 ' 

where the the symmetry with respect to the lattice mid point r/2 + 1 is enforced. In figure 2 we plot 
the pion effective mass squared at different time slices for quark masses am = 0.01 , 0.02, 0.03, 0.04, 0.05. 
Since effective masses show a flat behaviour one can pick the value at a given time slice. We have 
selected the values at the last time slice which display the largest errors. 
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Figure 2. Pion effective masses squared for different quark masses. 

In figure 3 we show the behaviour of the pion mass squared against the quark mass. The fuU 
line is the least squares fit of the data which gives the result: 

(am;,)^ = 0.007(17) + 3.5(4) (am) . 

The figure shows the extrapolated pion mass at zero quark mass, 0.007(17), which in dimensionful 
units gives a pion mass mj^ ~ 170(200) MeV. Here, we have assumed that, at this coupling, the 
inverse lattice spacing is ~ 2 GeV. At this accuracy, the pion mass is consistent to zero. 

To conclude, we have reviewed the recent efforts to revive the minimally doubled actions. 
We have made a preliminary calculation of the pion mass on a SU(3) background, which behaves 
as predicted by the chiral perturbation theory within the statistical error bars. The final results 
of this ongoing calculation will reveal any possible discrepancy to the chiral perturbation theory. 
One notices, however, that the invested computational effort to obtain these results is much smaller 
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than for the Ginsparg- Wilson fermions. Therefore, we conclude that the minimally doubled action 
presented here is worth exploring in the future. 
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Figure 3. Testing chiral perturbation theory with minimally doubled actions: pion mass squared against 
the bare quark mass. At zero quark mass, the extrapolated pion mass is ~ 170(200). 
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